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ABSTRACT 

' We select a sample of about 50,000 early- type galaxies from the Sloan Digital Sky 

Survey (SDSS), calibrate fitting formulae which correct for known problems with pho- 
' tometric reductions of extended objects, apply these corrections, and then measure a 

number of pairwise scaling relations in the corrected sample. We show that, because 
they are not seeing corrected, the use of Petrosian-based quantities in magnitude 
limited surveys leads to biases, and suggest that this is one reason why Petrosian- 
based analyses of BCGs have failed to find significant differences from the bulk of the 
Q ■ early-type population. These biases are not present when seeing-corrected parameters 

' derived from deVaucouleur fits are used. Most of the scaling relations we study show 

C/3 , evidence for curvature: the most luminous galaxies have smaller velocity dispersions, 

larger sizes, and fainter surface brightnesses than expected if there were no curva- 
ture. These statements remain true if we replace luminosities with stellar masses; they 
. suggest that dissipation is less important at the massive end. There is curvature in 

^ ' the dynamical to stellar mass relation as well: the ratio of dynamical to stellar mass 

CsJ I increases as stellar mass increases, but it curves upwards from this scaling both at 

04 ■ small and large stellar masses. In all cases, the curvature at low masses becomes ap- 

' parent when the sample becomes dominated by objects with stellar masses smaller 

, than 3 x lO^^M©. We quantify all these trends using second order polynomials; these 

■ generally provide significantly better description of the data than linear fits, except at 

the least luminous end. 
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1 INTRODUCTION Measurements of brightest cluster galaxy scaling rela- 

tions have shown them to be unusual (Malumuth & Kirshner 
Early-type galaxy observables, luminosities L, half-light ^gg^^ ^gg^. ^ ^ ^^^^^^^ ^gg^. p^^^^^^ ^ ^auer 1995). 

radu its, mean surface brightnesses Je, colors, and veloc- a,..-,-ii ■ -n , , , ,■ r 

And statistically sigmncant detections ot curvature m many 

scaling relations across the entire population have now been 



ity dispersions a, are strongly correlated with one another. 

These scaling relations are usually described as single power- j /ry i , i nnnr. t ^ i nnn^ r, j- ^ i 

^-i/K _i/4 n 4 7-. T-n s\ made (Zaritsky et al. 2006; Lauer et al. 2007; Bernardi et al. 

laws [Re OcL',aOcL',g — r0Ca-,Re0Cle ' ), SUg- „„„„ „ u * l onn-r t • -t l onno u 4- J 

^ , ' ^ , , , 2007; Desroches et al. 2007; Lm et al. 2008; but see von der 

gestmg a single formation mechanism across the population. t • i . i r>^or,\ mi ■ i <■ . i • • , 

° fr 1111 Linden et al. 2007). the mam goal oi this paper is to ex- 

However, galaxy formation models suggest that the bright- i ii i i • -j j u ii oi t^- i oi 

, , 1 r /T-N ploit the large sample size provided by the bloan Digital bky 

est galaxies m clusters had unusual formation histories (De „ /i ot^oon j_ i • • j_ r 

, Survey (hereafter bDbs) to make precision measurements ot 
Lucia et al. 2006; Almeida et al. 2008), so they may follow , . ,. , 

/T, 1 1 1 . 1 n'^^r 1 the curvature m these scaling relations, 

dinerent scaling laws (Boylan-Kolchm et al. 2005; Robert- 

son et al. 2006; Hopkins et al. 2008a). Formation histories. Section H describes how we select a sample of early-type 

and the importance of gaseous dissipation and/or gas rich 8^^^^^^^ ^^^^ ^^e SDSS. It also discusses the corrections we 

1 ijiu u j-CT ij J- apply to account for the fact that the SDSS photometric 

mergers are also expected to have been dinerent depending ^ 

on the mass range of the galaxy (e.g. Mihos & Hernquist reductions are unreliable for extended objects. These are 

1993; Naab et al. 2006; Hopkins et al. 2008b, Ciotti et al. particularly important since the curvature we would like to 

onn'7^ o • U4- ui j- i J -t measure is small (else it would have been seen in smaller 

2007). bo, one might reasonably expect to see departures ^ 

from the single power-law scahng relations, especially at the samples) , so the photometric and spectroscopic parameters 
extremes of the population. extremes of the sample must be reliable. 

Section [3] quantifies the curvature in a number of pair- 
wise scaling relations for this sample. It also shows that, in 

* E-mail: jhyde, bernardm@physics.upenn.edu the relations which involve luminosity, the curvature per- 
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sists if one replaces luminosities with stellar masses. A final 
section summarizes our findings. 

An Appendix discusses why, because they are not 
seeing-corrected, Petrosian based quantities are ill-suited for 
precision measurements in deep, magnitude-limited, ground- 
based datasets. It also shows why the Petrosian-based anal- 
ysis of BCGs by von der Linden et al. (2007) yielded anoma- 
lous results. 



2 THE SAMPLE 

We start from a sample of 376471 galaxies based on the 
Fourth Data Release (DR4) of the SDSS but with parame- 
ters updated to the SDSS DR6 values (Adelman-McCarthy 
et al. 2008). From this sample we extract 46410 early-type 
galaxies following the procedure described below. We use 
SDSS deVaucouleur magnitudes and sizes, model colors, 
and aperture corrected velocity dispersions to re/8 unless 
stated otherwise. The SDSS also outputs Petrosian magni- 
tudes and sizes. However, these are not seeing corrected, 
and Appendix A shows that this introduces systematic bi- 
ases, so we do not use them in what follows. Throughout, 
angular diameter and luminosity distances were computed 
from the measured redshifts assuming a Hubble constant of 
Hq = 70 km s~^ Mpc~^ in a geometrically flat background 
model dominated by a cosmological constant at the present 
time: (^o, Ao) = (0.3,0.7). 

About notation: we use to specify radius in kpc, and 
re to specify angular size in arcseconds. For surface bright- 
nesses, we use the following definition. 



/ie = -2.51ogio(Je) = -2.51ogio 




= m-f51ogio(re) -h2.51ogi(,(27r) - 10 logu,(l + z)(l) 

where m is the evolution, reddening, and k-corrected appar- 
ent magnitude. We use to denote the absolute magnitude 
in band x, and h'h to denote stellar mass in solar units. 

2.1 Selecting early- types 

To obtain a sample of early-type galaxies we first select 
the subset of galaxies which are very well-fit by a deVau- 
couleur profile in the g and r bands (p-band fracDev = 1 
and r-band fracDev = 1); this gives about 100603 objects. 
To avoid contamination by later-type galaxies we also re- 
quire the spectrum to be of "early-type", by setting eClass 
< (see SDSS documentation for details of this classifica- 
tion). This slightly reduces the sample, to 94934 galaxies. 
Since the SDSS spectroscopic survey is magnitude limited, 
we require that r— band deVaucouleur magnitudes satisfy 
14.5 < rrir < 17.5. (Spectra are actually taken for objects 
having Petrosian magnitude mr,pet < 17.77; our more con- 
servative cut is designed to account for the fact that the 
Petrosian quantity systematically underestimates the total 
light in a deVaucouleur profile.) This cut leaves 70417 galax- 
ies. Of these, 64492 have stellar mass estimates from Gallazzi 
et al. (2005). 

We would also like to study the velocity dispersions of 
these objects. One of the important differences between the 
SDSS-DR6 and previous releases is that the low velocity 
dispersions {a < 150 km s~^) were biased high; this has 




Log,o c'iDLspec2d [km s '] 

Figure 1. Comparison of the velocity dispersion measurements 
from the SDSS-DR6 and the IDLspec2d reduction. Vertical 
dashed line shows our cut at small cr ~ 60 km s~^. Symbols with 
error bars show the median value and its uncertainty, dashed and 
dotted curves show the regions which enclose 68% and 95% of the 
objects. 



been corrected in the DR6 release (see DR6 documentation, 
or discussion in Bernardi 2007). We compared the values 
given by the official SDSS-DR6 database and those com- 
puted in the Princeton reduction (lDLspec2d). The two sets 
of measurements still show weak systematic trends. The up- 
per panel of Figure [T] shows that at fixed o-dr6 the median 
values of criDLspec2d agree quite well with the aum values, 
except at large 0"dr6 > 320 km s~^, where (tdr6 is slightly 
larger. However, the bottom panel shows that, when com- 
pared at fixed criDLspGc2d, a systematic trend is more evident 
- especially at small CTiDLspcc2d < 120 km s~^. To minimize 
systematics, we decided to average the DR6 and IDLspec2d 
velocity dispersion measurements. 

At the low end, we select galaxies with cr > 60 km s~^ 
(see Section [3]4l Bernardi et al. 2003c and Hyde & Bernardi 
2008 for discussion of biases introduced by eliminating ob- 
jects based on their a). At the high end we select galax- 
ies with a < 400 km s~^ to avoid contamination from 
double/multiple superpositions (see Bernardi et al. 2006b, 
2008). The maximum cr of a single galaxy given in Bernardi 
et al. 2008, and confirmed by Salviander et al. 2008 using 
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Figure 2. Distribution of b/a in the r-band as a function of angular (left) and physical size (right). Contours show regions of equal 
probability density. Starting at maximum density, each line represents a factor of \/2 decrease in density. The change from solid to dotted 
lines marks the point which encloses 68% of the sample. There appears to be a separate population of small b/a < 0.6 objects at physical 
sizes smaller than a = I0h~^ kpc. 
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Figure 3. Histogram showing distribution of b/a in the g— (solid) 
and r— (dashed) bands. 

high resolution spectra, is about 430 km s~^. In addition, 
the SDSS-DR6 only reports velocity dispersions if the S/N 
in the spectrum in the restframe 4000 — 5700 A is larger 
than 10 or with the status flag equal to 4 (i.e. this tends to 
exclude galaxies with emission lines). To avoid introducing 
a bias from these cuts, we have also estimated velocity dis- 
persions for all the remaining objects. It turns out that this 
cut does not change the correlations studied here (nor those 
in the Fundamental Plane study of Hyde & Bernardi 2008) . 
The net result is to change our sample from 64492 to 64343 
objects. About 25 objects have colors which lie outside the 
range 0.4 < g — r < 1.3 which we do not believe, so we 
exclude them from the sample. 

Figure [5] suggests that we should make one more cut. 
The left hand panel shows the distribution of axis ratios 
6/a in the sample, as a function of the angular half light 
radius (left) and physical scale (right). There is clear evi- 
dence for two populations, particularly when plotted as a 
function of physical scale: one with b/a > 0.6 and the other 
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Figure 4. Dependence of axis-ratio b/a on a in the r— band for 
galaxies with b/a > 0.6 in the g— band. Contours show the joint 
distribution of b/a and a with the same conventions for spacing 
and line styles as in Figure |2] 



with b/a < 0.6. The objects with small b/a make up about 
twenty percent of our full early-type sample, but they are a 
larger fraction of the smaller fainter galaxies than they are of 
the bright. There are good physical reasons to suspect that 
these objects are a different population (e.g. rotational sup- 
port is necessary if b/a < 0.6), so, in what follows, we remove 
all objects with b/a < 0.6 in the g—hand from the sample. 
This leaves 51379 objects. Many of the figures which follow 
show quantities in the r— band, for which this cut appears 
less sharp. Figure [3] shows that the distribution of b/a in the 
two bands is very similar: the small differences between b/a 
in g and r is due, in part, to measurement error. 

Figure |4] shows b/a vs a in the r— band in our sample 
after applying this cut. Comparison with Figure [5] shows 
that the 'second' population at low b/a has been removed. 
In the sample which remains, there is a weak tendency for 
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the largest galaxies to have slightly smaller b/a. The impli- 
cations are discussed further in Bernardi et al. (2008) who 
show that the mean h/a drops slightly at the highest lumi- 
nosities and a. 

Finally, we must account for the fact that objects were 
brighter in the past by 0.9 z mags because of stellar evolution 
(e.g. Bernardi et al. 2003a). This, with fc-corrections, adjusts 
slightly the actual values of the magnitude limit which we 
should use when computing effective volumes. In fact, these 
effects work in opposite directions, so the net effect is small. 
Ifowever, since our goal is to quantify small effects in a large 
sample, it is necessary to do this carefully. The net result is 
to reduce the sample size by about ten percent, to 46410. 
As a check that this has been done correctly, we perform 
the test suggested by Schmidt (1968). If Vi is the survey 
volume between object i and the observer, and Vmax.i is the 
total survey volume over which the object could have been 
seen, then the average value of (V/Knax) should be 0.5: we 
find {V /V-mia) = 0.499. (The luminosity evolution is slightly 
larger than, but statistically consistent with, the 0.852 mags 
reported by Bernardi et al. 2003b.) 

These cuts are similar to those described by Bernardi 
et al. (2003a, 2006a), who provide further details about the 
motivation for each cut, except for: i) the cut on h/a; ii) 
the inclusion of velocity dispersions from low S/N spectra 
or with the status flag not-equal to 4; and iii) the inclu- 
sion of velocity dispersions at 60 < cr < 90 km s~^. In 
addition, for the present study, we have chosen to be more 
conservative. Previously, we required Petrosian magnitudes 
14.5 < mr.Pot < 17.75. Changing to a brighter limit makes 
our final sample size considerably smaller than if we had 
used the Bernardi et al. (2006a) selection. In addition, we 
previously required f racDev > 0.8 in the r— band; we now 
require fracDev = 1 in r as well as in g, because non-early- 
type features are expected to be more obvious in the g— 
band. This more conservative choice for f racDev eliminated 
about 20000 objects (doing the selection based on g but not 
r makes little difference, because requiring f racDev = 1 is 
quite stringent). To quantify the effect of these additional 
cuts, we have applied them to the Bernardi et al. (2006a) 
sample, and found that the comoving number density of ob- 
jects is reduced to about 0.4 times that in Bernardi et al. 
(2006a). None of the results which follow are sensitive to the 
value of the comoving number density. 



2.2 Corrections to photometry 

We must address another issue before proceeding. This is be- 
cause the SDSS reductions are known to suffer from sky sub- 
traction problems, particularly for large objects (Adelman- 
McCarthy et al. 2008). To illustrate. Figure [5] compares 
SDSS photometric reductions with those of our own code, 
GALMORPH, for a subset of the fuU sample (~ 5500 galax- 
ies used by Bernardi et al. 2003a plus ~ 180 Bright Clus- 
ter Galaxies analyzed by Bernardi et al. 2007). The GAL- 
MORPH reductions do not suffer from the sky subtraction 
problem. This figure shows that while the two pipelines are 
in excellent agreement for small objects, the SDSS underes- 
timates the sizes and brightnesses of large objects. However, 
the quantity IP = log lO(re) — 0.3 /ie, identified by Saglia et 
al. (2001) (they called it FP), is not substantially changed. 
In all the panels, symbols with error bars show the mean 
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Figure 5. Comparison of GALMORPH and SDSS photometric 
reductions. Although the two are in good agreement at small re, 
where the bulk of the objects lie, the SDSS underestimates the 
total flux (top) and the half-light radius (middle) of large objects. 
However, both pipelines return consistent values of the quantity 
logio(''e)-0.3 ^le which Saglia et al. (2001) call FP. Symbols with 
error bars show the mean relation and its uncertainty, dashed 
and dotted curves show the regions which enclose 68% and 95% 
of the objects, and smooth solid curves show the fits given by 
equations (|4]l and JSj. 
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[arcsec] 

Figure 6. Distribution of sizes obtained by the SDSS photometric 
reduction (sohd) and corrected following Eq. [5] (dotted). Although 
the two are in good agreement at small rg, where the bulk of the 
objects lie, the corrected sizes are larger for the small fraction of 
large . 



relation and the error on the mean, dashed curves show the 
region encompassing 68% of the objects, and dotted curves 
enclose 95%. 

Unfortunately, it is computationally expensive to run 
GALMORPH on the entire DR6 data release. Therefore, we 
have fit smooth curves to the trends shown in Figure [5] (solid 
curves) and we use these fits to correct the SDSS reductions 
as follows. Given rsoss we set 



re = rsDSS + Ar-fit(rsDSs) 
m = msDSS + Amflt(rsDSs), 



(2) 
(3) 



where rsDSS = rjevy - is the SDSS-measured deVau- 

couleur radius expressed as a geometric mean (-\/a6) of the 
semimajor (a) and semiminor (6) axes of a half-light con- 
taining ellipse, rdev, the semimajor axis length, and -, 
the axis ratio, are obtained from the SDSS catalogs where 
they are referred to as "deVRAD" and "devAB", respec- 
tively. Arflt(T'SDSs) is zero if rsoss < 2 arcseconds and 
Amfit(rsDSs) is zero if tsdss < 1.5 arcseconds, otherwise 

Amfit(rsDSs) 



mags 
Arfit(rsDss) 



= 0.024279 



= 0.181571 



71.1734 ^ 

rSDSS _|_ 



26.5 

/ rsDss y 
V 3.9165/ 



(4) 
(5) 



4.5213 V 3.9165/ 
We propagate the errors similarly. Note that 
—0.6 In(lO) rsDSS Amflt provides a good approxima- 
tion to Arflt, as one might expect given the bottom panel in 
Figure [H it slightly underestimates Arat at large rsDSS- We 
have tested these corrections, and found them applicable to 
the SDSS g, r, and i bands. Throughout the paper we will 
denote angular size in arcseconds with and physical size 
in h~^ kpc with Re. These values refer to geometric-mean 
deVaucouleur radii, corrected as described in this section 
and Section [2.31 

Later in this paper, we will use stellar mass estimates 
from Gallazzi et al. (2005). These are derived by estimating 
a stellar mass-to-light ratio, and then multiplying by the 
estimated luminosity. Since our corrections to the appar- 
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Figure 7. Observed (top) and corrected (bottom) sizes in the 
SDSS derived from deVaucouleur fits to the light profiles shown 
as a function of r-band luminosity Mr- (In all cases, A logj^Q Re = 
logj^Q Re + 0.22 Mr + 4.24.) g-, r-, i-, and z-bands are shown from 
top to bottom in the upper panel. The bottom panel omits the 
2-band. In each panel, for each band, different lines show data 
from the following redshift bins: 0.07 < z < 0.1, 0.1 < z < 0.13, 
0.15 < z < 0.18, 0.22 < z < 0.25 and 0.25 < z < 0.35. 



ent magnitudes will affect these luminosities (they increase 
slightly on average), we add equation (3) to Gallazzi et al. 
stellar mass estimates (accounting for the conversion from 
luminosity to magnitude, m = — 2.51ogjg L). The net effect 
is to slightly increase some of the stellar masses, but we have 
not included a plot showing this increase. 

Figure |S] shows the distribution of the effective radii 
of the 46410 early-type galaxies (i.e., whatever their b/a 
value) before and after applying the correction given in equa- 
tion (O. Although the two are in good agreement at small 
re, where the bulk of the objects lie, the corrected sizes are 
larger for the small fraction of objects with large Ve- 

2.3 Corrections to sizes 

We make one final correction to the sizes, which is aimed 
at accounting for the fact that the early-type galaxies have 
color gradients: on average, their optical half-light radii are 
larger in bluer bands. If not accounted for, a population 
of intrinsically identical objects will appear to be slightly 
but systematically larger if they are at higher redshift. Our 



6 J. Hyde & M. Bernardi 



sample is large enough that we must correct for this effect, 
the moral equivalent of the fc-correction to the luminosities. 
We do so following Bernardi et al. (2003a). 

We estimate the rest-frame radius in each band by inter- 
polating the observed radii in adjacent bands. For example, 
to estimate the rest-frame g-band radius, we use the follow- 
ing expression: 



R. 



_ (1 + Z)\g - \r I 
■e,g,rest — , , \J^e.g ,ohs 



R 



e.r.obs 



(6) 

where \g,r,i,z = {4686A, 6165A, 748lA, 893lA} are the 
average wavelengths of the SDSS filters, and z is the 
spectroscopically-determined redshift. 

Figure [7] shows why this correction is necessary. It shows 
the sizes of objects in different bands as a function of r- 
band {k- and evolution-corrected) luminosity. To reduce the 
range of sizes, we have subtracted-off a luminosity depen- 
dent factor: we actually show Alogj^Q Re = logio(^e/kpc) -I- 
0.221 Air + 4.239 (the reason for the exact choice of param- 
eters will become clear shortly - see Table 1 - but note 
that, for the present purpose, the exact choice is not im- 
portant). The upper panel shows the observed sizes: from 
top to bottom, the sets of symbols and lines show g-, r-, i-, 
and z-bands. For each band, different lines show data from 
a number of redshift bins: 0.07 < z < 0.1, 0.1 < 2 < 0.13, 
0.15 < z < 0.18, 0.22 < z < 0.25 and 0.25 < z < 0.35. In 
any redshift bin, the sizes are clearly larger in g- and smaller 
in i- than they are in the r-band. The bottom panel shows 
rest-frame sizes, for g-, r-, and i-bands. We omit the z-band 
since longer wavelength observations would be necessary to 
reconstruct the z-band rest-frame size. The rest-frame sizes 
are indeed larger in the bluer bands, with the difference per- 
haps slightly smaller at large luminosity. The redshift depen- 
dence of the size-luminosity relation is not apparent when 
using the the observed radii (upper-panel), but can be seen 
when using the rest-frame radii (lower panel); this is studied 
further in Bernardi (2009). 

The SDSS also outputs non-parametric Petrosian sizes. 
However, in contrast to the deVaucouleur-fitted quantities, 
these sizes are not corrected for the effects of seeing. The 
Appendix shows that, in contrast to the deVaucouleur sizes, 
the Petrosian sizes of objects at higher redshift are systemat- 
ically larger than those from the model based fits. This is not 
surprising if seeing has compromised the Petrosian-based 
measurements, suggesting that, if this is not accounted for, 
then the use of Petrosian sizes limits or biases the precision 
measurements which large sample sizes would otherwise al- 
low. 



3 CURVATURE IN SCALING RELATIONS 

We now turn to measurements of a number of scaling rela- 
tions. It turns out that curvature is often even more obvious 
if we replace luminosity with stellar mass, so we will often 
show such relations side by side. Unless we specify otherwise, 
the luminosity is always from the r-band. We will some- 
times use a shorthand for the r-band quantities: R, V, I, M 
for logio Re, logio a, fie, Mr. 
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Figure 10. As for previous figure, but now for the ratio of dy- 
namical mass to light versus luminosity (top) and dynamical to 
stellar mass versus A/» (middle) and M^y^/L or M^y^/M, versus 
cr (bottom). 
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Table 1. Coefficients of best-fitting relations of the form |X) = po + Pi + P2 to pairwise scaling relations. Fits were made to 
the binned points, not to the objects themselves. Linear fits (p2 = 0) were made to the galaxies and restricted to the range 10.5 < 
l°glo(^*/^0) ^ ^^-5 and —23 < < —20.5. The errors on the fitted coefficients are random errors: they depend on slope of the 
relation, its scatter, and the sample size. These are smaller than those produced by systematic effects (e.g., using a"DR6 o'' flDLspoc2d 
rather than their average). Typical systematics errors are a few times larger than the random errors. 
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Figure 8. Half-light radius vs stellar mass and luminosity (left and right). Symbols with error bars show the median in small mass or 
luminosity bins, and dashed and dotted lines show the regions which contain 68% and 95% of the objects in each bin. Curves show fits of 
the form (Y|X) = po +Pl ^ + P2 to these relations; best-fit coefficients are provided in Table [T] Curved fits were made to the binned 
counts (symbols with error bars), rather than to the objects themselves. Straight lines show linear fits to these relations; these were made 
to the galaxies and restricted to the range 10.5 < \og-^Q(Mf / Mq) < 11.5 and —23 < Mr < —20.5. Bottom panels show residuals from 
the Vmix-weighted linear (connected thin line) and quadratic fits (filled circles). 



8 J. Hyde & M. Bernardi 




3.1 Curvature in pairwise relations 

A simple test for curvature is as follows. In a magnitude 
limited survey, the more luminous objects are seen out to 
bigger volumes than the fainter objects, meaning that the 
ratio of the number of luminous to faint objects seen in such 
a survey is larger than the ratio of the true number densities 
of these objects. One can account for this over-abundance by 
weighting each object by l/ymax(i), the inverse of the vol- 
ume over which it could have been seen; this down-weights 
luminous galaxies relative to fainter ones. If a scaling rela- 
tion is a straight line, e.g., if the mean log(size) increases 
linearly as log(luminosity) increases, then the slope of the 
regression of log R on log L will not depend on whether or 
not one includes this weighting term. (This will also be true 
for log a on A/, etc.) However, if the intrinsic scaling re- 
lation is curved, then the slope of the regression line will 
depend on which weighting was used. This happens to be 
true in our dataset, indicating that the relations are curved: 
the slope of the log R - M relation is -0.221 ± 0.001 when 
weighting by KiTaxi but —0.241 when not. These numbers 
are -0.119 ± 0.001 and -0.117 for the log a - M relation. 
The difference is more dramatic for the /le — M relation: 
-0.104 ± 0.004 and -0.207. Here, we reported the uncer- 
tainties due to random errors; uncertainties in the param- 
eters due to systematics errors are a few times larger (see 
Section 

3.2 The size-luminosity relation 

Recent work shows that the correlation between size and 
luminosity relation has evolved significantly since z = 2 
(e.g. Cimatti et al. 2008; van Dokkum et al. 2008). Bernardi 
(2009) shows that the sizes of luminous early-type galaxies 
(Mr ~ —23 mag) are still evolving at low redshift and that 
satellites are on average ~ 8% smaller than central galax- 
ies. However, this last result is not seen for lower luminosity 
early- type galaxies (e.g. Weinmann et al. 2009). Due to this 
recent interest, we study this relation first, before consider- 
ing others. 



The symbols in the panels on the right of Figures|8]and[9] 
show the curvature in the size-luminosity and cr-luminosity 
relations, respectively (e.g. Lauer et al. 2007; Bernardi et 
al. 2007). This curvature is also evident when shown as a 
function of stellar mass (panels on the left), so curvature in 
the Mt/L relation is not the only cause. 

The format of this figure is similar to many that follow. 
The filled circles in the top panel show the median size in 
a number of narrow bins in luminosity, when the objects 
are weighted by V^nTax- Dashed and dotted curves show the 
regions which include 68% and 95% of the weighted counts. 
For statistics at fixed luminosity (as in this case), the Vmax 
weighting makes no difference because, at fixed luminosity, 
all galaxies have the same weight. In Figure [12] which shows 
the R — a and R — jie relations, the difference is significant. 
There, we use open squares to show the median weighted 
count when no Vmax weighting is used. 

Straight solid and dashed lines show the result of fitting 
straight lines to the data, with and without Vmax weighting. 
For distributions at fixed luminosity (as in this case), the 
two should be the same if the underlying scaling relation is 
not curved. If it is, then the fit to the unweighted points will 
refiect the slope of the relation at higher luminosities. In the 
case of the size- luminosity relation (Figure 8), the dashed 
line is steeper than the solid, consistent with the steepening 
of the relation at large L shown by the symbols. In this case 
only, we also show the linear fits to the R — L and R — M* 
relations reported by Shen et al. (2003) (see their Table 1, 
Fig. 6 and Table 1, Fig. 11 respectively). This shows that 
their fits are close to ours when we ignore the Knax weight- 
ing, but note that their sample is selected rather differently, 
and their sizes are from Sersic, rather than deVaucouleur 
fits to the light profile. 

To quantify the curvature, we fit 2nd order polynomials 
to these and a number of other scaling relations (which fol- 
low). In all cases our fits are slightly non-standard because, 
to emphasize the curvature in these relations, we would like 
the fits to be sensitive to the tails of the distribution. There- 
fore, we have fit to the binned counts (i.e., the symbols) 
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shown in the Figures, rather than to the objects themselves. 
For fits at fixed L, this is equivalent to weighting each object, 
not just by V-^L), but by [Knax(L)0(L)]-' = [Nobs{L)]-\ 
In efltect, this upweights the tails of the distribution. The 
fitting minimizes x^i defined as the sum of the squared dis- 
tances to the binned points shown. If, instead, we weighted 
each of the binned points by (the inverse of its) error bar 
when defining x^, then, because this additional weight is 
proportional to the number of counts in the bin, this would 
be the equivalent of weighting each object, rather than each 
binned point, equally. 

The bottom panel shows residuals from the linear (V^^^ 
weighted) and quadratic fits; in all cases, the residuals from 
the quadratic fits show fewer, if any, trends. Table 1 summa- 
rizes the results of these fits. Although the coefficients of the 
quadratic term appear to be very different from those when 
simply fitting a line, this is because we are reporting fits of 
the form = po + pi X + p2 . Had we removed the 

mean values first, and fit {y\x) = po + Pi x + p2 [x^ — {x^)] 
instead (where x = X — (X), etc.) then pi would be very 
similar to linear fit, indicating that the curvature is small. 

To quantify if a quadratic is a significantly better fit 
than a straight line, we fit both to the binned counts. We 
then compare Xquad/(^bms-3) with Xune/(^bins-2), where 
Xquad and Xune dcuote the minimized values of x^, and the 
fits were to A'^bins binned points. In all cases, this reduced 
value of for tli6 quadratic is much closer to unity than 
for the linear fits. (Note that we do not use the coefficients 
of the linear fits to the unbinned counts that are shown in 
Table 1, since, for this test, we want to treat both the linear 
and quadratic fits equally. Of course, using these coefficients 
when computing Xuno does not change our conclusions, since 
they can, at best, produce the same value of Xunc which we 
have just described.) 

Finally, a word on the errors on the fitted coefficients is 
in order. The numbers we quote are random errors: they de- 
pend on slope of the relation, its scatter, and the sample size. 
These are smaller than systematic effects: e.g., using (Jdr6 
or (TiDLspcc2d rather than their average makes a difference 
which is larger than the random error. Typical systematics 
errors are a few times larger than random errors - therefore, 
it is important to separate the two types of error. 

Figure |8] shows that, compared to a single power law, 
the R — L relation curves significantly towards larger sizes 
at large L or M«, consistent with previous work. At small 
L or M. {Mr < -19.5 or log-^^o M,/Mq < 10.2), the data 
appear to scatter slightly downwards from the quadratic fit. 
Although there are few objects in this tail, so the measure- 
ment is noisier, it is possible that this indicates that the 
sample is slightly contaminated at small L. We will return 
to this shortly. 

Before we move on to other scaling relations, we note 
that results based on Petrosian quantities are shown in the 
Appendix, where we also discuss why we do not consider 
them further. 

3.3 Other scaling relations 

Figure |9] shows the a — L relation in this sample. This re- 
lation is actually rather well described by a single power- 
law, except at A'lr < —23 (where the mean value of 
logjo((T/km s~^) > 2.4) where it curves slightly downwards. 



The flattening of the a — L relation is consistent with previ- 
ous work, but notice that this effect is even more pronounced 
for the a — A4, relation. 

Comparison with Figure |S] shows that, at large M» and 
L, the R — L and a — L relations curve in opposite senses. 
Since dynamical mass Mdyn oc _Ra^, it is interesting to check 
if the Mdyn — L relation is well-fit by a simple power-law. 
We do this in Figure [TOl where we set Mdyn = bR^a^ /G — 
4.65 X lO^°/i-^M0 (7?e//i"^kpc) (o-/200 km s'^f. The top 
panel shows that, at both small and large L, the Mdyn — L 
relation curves upwards from the Ra^ oc L"'^ scaling (solid 
line). The middle panel shows that curvature remains if 
one replaces L with M*, indicating that more than stellar 
population related effects are responsible. The slight rise at 
small masses is not implausible: star formation is expected 
to be less efficient at small masses. However, we view this 
with caution: the velocity dispersions at the small-mass end 
are more uncertain (e.g. Bernardi 2007). There is an aver- 
age trend for Mdyn/M« to increase with M«, although it is 
weak: {Mciyr,/M^\Md.y„) oc M° °«2±°-°°'' (the error 0.06 on 
the slope was computed accounting for systematics errors 

- the uncertainty from random errors is smaller ~ 0.002). 
The correlation between Mdyn/M» with Mdyn is stronger: 
{Miy^/M,\Miyr,) OC M°^ij±0'°i (see Figure 11 in Hyde & 
Bernardi 2008). If there were no scatter around this rela- 
tion, then we would expect {Mdyn/Mt\Mt) oc O-i'^/^ss ^ 
M^'^; because there is scatter, this scaling is shallower, 
oc jvf0 06±o.oi gygjj jf ^jjjg relation were a simple power 
law without curvature, it would indicate that stars make 
up a smaller fraction of the total mass of a galaxy at large 
masses. This provides an important piece of information for 
adiabatic contraction based models of scaling laws (e.g. Pad- 
manabhan et al. 2004; Lintott et al. 2005). 

To connect with previous work, the bottom panel shows 
how Afdyn/i (triangles) and Mdyn/M« (filled circles) scale 
with cr; both relations are slightly curved. Except at the 
largest a, our data are relatively well described by the curved 
MdynZ-f/ — o" relation reported by Zaritsky et al. (2006) (we 
have shifted our measurements downwards by 0.44 dex be- 
cause our data are in r whereas their fit was in /, i.e. we use 
r — I — 1.1. Note that we have also subtracted —0.3 from 
their fit since they used the effective light Le = L/2, while 
we use the total light L). The fact that Mdyu/L is a steeper 
function of a than is Mdyn/M* can be understood as follows. 
First, note that Mdyn/L = (Mdyn/M*) {M,/L). Then, note 
that Mf/L increases with increasing g — r color (e.g. Bell et 
al. 2004). However, g — r color is strongly correlated with a 

- large a implies redder colors (Bernardi et al. 2005), and so 
Mt/L increases with a. Therefore, Mdyn/L increases with a 
because Mdyn/M* does and because M,/L does so as well. 

As a final study of curvature in relations which involve 
luminosity, we now turn to the fJ,e — L and /i* — M« relations. 
(We define fj., = -2.5 logio(M./MQ) + 5 logio(i?e/kpc) + 
2.5 log]^Q(27r); this is the stellar mass surface brightness 
within the half light radius.) Figure [TT] shows that in this 
case too, there is signiflcant curvature. However, the panel 
on the left shows that at M* < 3 x 10^° Mq, the /x* - M, 
relation becomes rather well flt by the linear relation, al- 
though it also becomes significantly noisier. A look back at 
the other relations which involve M* shows that they too 
become less well-defined at small M*. This happens to be 



10 J. Hyde & M. Bernardi 




9.5 10.0 10.5 11.0 11.5 12.0 12.5 -18 -19 -20 -21 -22 -23 -24 

Log,o M. [Ms„„] M, [mag] 



Figure 11. As for previous figure, but for stellar mass surface- brightness vs stellar mass (left) and surface-brightness vs luminosity 
(right). 
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Figure 12. As for previous figure, but now for the size-velocity dispersion and size-surface-brightness relations. 



the same mass scale which Kauffmann et al. (2003) identify 
as being special. 

It is possible that our early-type sample is contaminated 
by a different population at the low mass end - despite the 
fact that we have already tried to reduce such an effect by 
removing objects with b/a < 0.6 and selecting galaxies with 
the g and r-band fracDev = 1. If we include those objects 
with g and r-band fracDev > 0.8 and/or do not remove 
objects with small b/a, then the quadratic remains a good 
fit even at small M, or L. This is also true if we simply use 
the cuts given by Kauffmann et al. (2003): R90/R50 > 2.86, 
R50 > 1.6 and fi^o < 23 mag/arcsec^ in the r-band. We 
will return to this in the next subsection, but note that the 
curvature at the luminous, massive end of the sample, is 
highly significant. 

We turn now to the scaling relations which play a fun- 



damental role in the Fundamental Plane: the Re ~ a and 
Re — 1-i.e relations. Figure \n\ shows that both these relations 
are curved. Notice that accounting for selection effects is 
important - the filled and empty symbols, which include or 
ignore the V^^x weight, trace very different relations. To first 
order, the zero-points of the two relations are more strongly 
affected than is the slope. Since it is the zero-point of the 
FP which is used to estimate evolution in small high-redshift 
samples. Figure [121 suggests that, without due care, one may 
simply be measuring selection effects (a point also made by 
Bernardi et al. 2003c). 

Neither of these relations is as well-fit by a quadratic at 
the extremes; e.g., when weighted by Vmax, the Re — He rela- 
tion curves away significantly from the quadratic at both 
large and small /Je (see also Nigoche-Netro et al. 2008). 
Curvature at large fie (or a) is not surprising: we already 
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Figure 13. PCA eClass (top) and angular size (bottom) as a 
function of stellar mass. Later- type galaxies have less negative 
eClass values. 



know that BCGs, and more generally, large mass galaxies 
(Bernardi et al. 2008) , follow different scaling relations. The 
flattening at small a is perhaps more surprising; this is in 
the regime where the SDSS velocity dispersions are most 
suspect, so one might worry that some of the effect is from 
measurement errors scattering objects to smaller a. How- 
ever, there is a similar flattening of the Re — fie relation at 
/ie < 18. What causes this? 



3.4 Curvature from contamination? 

In the previous subsection we raised the possibility that the 
sample may be contaminated at low M, . With this in mind, 
we performed a visual inspection of the objects with small 
M*; this did not reveal any peculiarities. However, the top 
panel of Figure [TJ] shows the results of a PCA analysis of 
the spectra. A small departure from early- towards later- 
type values is seen at small- M«, but note that true late- 
types have eClass values which are greater than zero. If we 
assume that some fraction / of the sample has eClass — 
0, and the rest has eClass = —0.13, then to find eClass — 
—0.11 requires / « 0.15. The true contamination is likely to 
be smaller (since eClass > for late-types). 

Alternatively, one might have worried that this is an 
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Figure 14. Similar to previous figures, but for the color- 
magnitude and color-stellar mass relations. 



aperture effect associated with abundance gradients - the 
spectra are from fibers which take light from a fixed angular 
radius, so, for smaller galaxies, the light from the inner bulge 
contributes a decreasing fraction of the total light in the 
fiber. However, a plot of Ve vs M* is approximately constant 
at small M» (bottom panel of Figure [13)) . suggesting that 
aperture effects are not to blame. 

Finally, we note that the color-stellar mass and color- 
magnitude relations also show a small curvature, towards 
bluer colors, at small L or M, . The curvature in these re- 
lations is enhanced if we allow the full range of b/a in our 
sample (rather than removing objects with b/a < 0.6) and 
if objects with g and r-band fracDev > 0.8 are included 
(rather than selecting galaxies with fracDev = 1). This, 
further circumstantial evidence for a small amount of con- 
tamination, is shown in Figure [14] 

We have also studied what happens if we remove objects 
with a < 90 km s~^ . This is prompted by Figure [12] and also 
by the fact that Bernardi et al. (2003a) excluded such ob- 
jects from their study of SDSS early-types (on the grounds 
that this is the regime in which the SDSS dispersions are 
suspect). Figure [15] shows that removing such objects has a 
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dramatic effect on the scaling relations which involve surface 
brightness, since this preferentially removes objects with 
large sizes for their luminosities (or stellar masses). As a re- 
sult, at the small mass end the size-luminosity (stellar mass) 
scaling relation leans towards smaller sizes. On the other 
hand, the a — L relation now flattens significantly at small 
L, because the small as at small L have been removed (there 
are essentially no objects with small a at large L), and the 
luminosity-/ie fiattens significantly. In addition, M^y^/M, 
vs M» becomes more curved because, for a given M«, one 
is removing small Mdyn- This leaves only large values of 
M^yn/M, sit small M*, making the relation curve upwards 
steeply. See Hyde & Bernardi (2008) for further discussion 
of the dramatic effects that cuts in a can have on the Fun- 
damental Plane (also see Bernardi et al. 2003c). 



4 DISCUSSION 

We used our own photometric reductions (GALMORPH) of 
about 6000 early-type galaxies from the SDSS to calibrate 
corrections to SDSS photometry (equations |3}[5]) which are 
most necessary for photometry of extended objects (Fig- 
ures [5] and |6]) . We applied these corrections to a larger sam- 
ple of about 50000 early-types, and then analyzed a number 
of galaxy scaling relations in the sample. Selection of the 
sample, which is described in some detail in Section |2] was 
more conservative than in previous work based on SDSS 
data. 

Small but statistically significant curvature was found 
for all the relations (Figures [8 HT2ll . Table [T] lists the coef- 
ficients of best-fitting second-order polynomials which pro- 
vide a concise way of describing this curvature. The Table 
also provides the coefficients of linear fits, for comparison. 
Whereas curvature at large luminosities/stellar masses is 
expected - BCGs are known to be a special population - 
some of the scaling relations in our sample show curvature 
at small masses as well. The critical mass scale is about 
M. « 3 X 10^° A/0. Whereas we see some evidence from the 
colors and spectra that objects below this mass scale are of 
slightly later type (Figures 1131 and I14p . a visual inspection 
of the images shows no obvious peculiarities. 

Our analysis also showed that the ratio of dynamical to 
stellar mass increases at large masses (Figure [TO)): this is a 
useful constraint on models of early-type galaxy formation. 
In addition, the R — a and i? — /ie relations were shown 
to be rather sensitive to selection effects (Figure I12p ; this 
matters for analysis of the Fundamental Plane. The ques- 
tion of whether or not the Fundamental Plane is curved 
or warped is addressed in a companion paper (Hyde & 
Bernardi 2008). Finally, we showed that seeing biases the 
scaling relations associated with Petrosian-based quantities 
(Appendix A), making them ill-suited for precision analyses 
in large ground-based datasets. 
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Figure 15. Effect on scaling relations of removing objects with cr < 90 km s~^. Smooth solid lines show the linear fits to the full sample 
(taken from previous figures); symbols show the scaling relations when objects with cr < 90 km have been removed; dashed and 
dotted lines enclose 68% and 95% of the objects. Parabolas show the quadratic fits to the symbols. 
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APPENDIX A: PETROSIAN-RELATED 
SCALING RELATIONS ARE BIASED BY 
SEEING 

The main text uses sizes and luminosities obtained from fit- 
ting deVaucouleur profiles to the images. We chose these 
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Figure Al. Residuals from tlie size-luminosity relation in the 
SDSS r-band, derived from deVaucouleur (top) and Petrosian 
(bottom) quantities; only the former are corrected for the effects 
of seeing. Different sets of symbols in each panel show this rela- 
tion at different redshifts; the relation is essentially independent 
of redshift in the top panel, but the high redshift objects appear 
to have larger sizes in the bottom panel. 

over the Petrosian quantities (R50 and petromag) because 
these latter are not corrected for the effects of seeing. To 
illustrate that seeing matters, Figure lAll shows the size- 
luminosity relation in six redshift bins (0 < z < 0.07, 
0.07 < z < 0.1, 0.1 < z < 0.13, 0.15 < z < 0.18, 
0.22 < z < 0.25 and 0.25 < z < 0.35, although the lu- 
minosities are k- and e-corrected to z = 0). The top panel 
shows deVaucouleur- and bottom shows Petrosian-based re- 
sults. In the top panel, the size-luminosity relations in the 
different redshift bins lie on top of one another. In the bot- 
tom panel, however, the high redshift relations are offset to 
larger sizes, suggesting that the Petrosian sizes have been 
biased high by seeing. 

Figure \K2\ presents other evidence that seeing matters. 
It shows the Petrosian concentration (the ratio of the ra- 
dius which contains 90% of the Petrosian hght, to that 
which contains 50%) as a function of redshift, for a few bins 
in (evolution-corrected) luminosity (—21 < Afr < —20.5, 
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Figure A2. Petrosian concentration as a function of distance 
from observer for a few bins in (evolution corrected) luminosity 
(-21 < Mr < -20.5, -22 < Mr < -21.5, -23 < Mr < -22.5, 
and —23.5 < Mr < —23) The higher redshift objects have smaller 
concentrations, because R50 is increasingly affected by seeing. 



-22 < Mr < -21.5, -23 < M,. < -22.5, and -23.5 < 
Air < —23). The catalog is apparent magnitude limited, so 
the most luminous bin extends to highest redshift. At fixed 
L, the higher redshift objects appear to be much less concen- 
trated, but this is not a physical effect: it happens because 
R50 is increasingly affected by seeing. 

In a magnitude limited survey such as the SDSS, the 
more luminous objects are seen preferentially at larger red- 
shifts. If not accounted for, seeing effects would bias the 
Petrosian size-luminosity relation (any scaling relations for 
which size or luminosity matter would also be affected) . This 
can lead to important differences in what one concludes from 
the data. 

For example, in their study of the R—L relation, von der 
Linden et al. (2007) used Petrosian quantities. They argued 
that BCGs traced essentially the same R — L relation as the 
bulk of the early-type galaxy population. This contradicted 
Lauer et al. (2007) and Bernardi et al. (2007) who found, on 
the basis of seeing-corrected fits, that the BCG relation was 
substantially steeper, and used this to draw important con- 
clusions about BCG formation histories. (Steeper relations 
for BCGs have since also been found by Liu et al. 2008). 

The bottom panel of Figure lAll shows that only when 
one stacks all redshift bins together (and ignores the obvi- 
ous offset from one redshift bin to the next) does one find 
a Petrosian R — L relation that is similar to that obtained 
from seeing-corrected fits; the relation in any one redshift 
bin is significantly shallower. However, at the highest red- 
shifts, where only the most luminous galaxies contribute, 
i.e., in the luminosity regime which is most likely to be dom- 
inated by BCGs, the R — L relation is indeed substantially 
steeper. Moreover, its slope becomes similar to that of the 
(incorrectly!) stacked sample, leading to the conclusion that 
BCGs have the same slope as the bulk of the population. In 
effect, using Petrosian quantities without appropriate care 
for the fact that they are not seeing-corrected, led von der 
Linden et al. to confusion. For this reason we have not used 
any Petrosian-based quantities in our analysis, and we cau- 
tion against their use in general. 



This is not to say that we think deVaucouleur-based 
fits are ideal. One might legitimately ask how the curva- 
ture we quantify in the main text changes if we fit to the 
more general Sersic profile. Fitting two-dimensional Sersic 
profiles to the SDSS photometry is well beyond the scope of 
this work. However, in principle, we could have done this ap- 
proximately as follows. Graham et al. (2005) provide a pre- 
scription for transforming from Petrosian to Sersic quanti- 
ties. So we could have used these transformations as a proxy 
for actual Sersic fits, and then studied the associated (PSer- 
sic?!) scaling relations. In fact, this was done by Desroches 
et al. (2007). We did not do this here because Figure lAll 
suggests that, because the Petrosian based quantities have 
not been seeing-corrected, this would yield biased results. 
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